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Abstract 

We investigate the vacuum polarization effects associated with a charged massless spin- 
1/2 field in a higher-dimensional space-time, induced by a composite topological defect. The 
defect is constituted by a global monopole living on a three-brane and two-dimensional con- 
ical space transverse to the latter. In addition, we assume the presence of an extra magnetic 
flux along the core of the conical space. The heat kernel and the Feynman Green function are 
presented in the form of a sum of two terms. The first one corresponds to the contribution 
coming from the bulk with global monopole in the absence of conical structure of the orthog- 
onal two-space, and the second one is induced by this structure and the magnetic flux. We 
explicitly evaluate the part in the vacuum expectation value of the energy-momentum tensor 
induced by the flux carrying conical structure. As in pure cosmic string geometries, only the 
fractional part of the ratio of the magnetic flux to flux quantum leads to non-trivial effects. 
The vacuum energy-momentum tensor is an even function of this parameter. We show that 
for strong gravitational fields corresponding to large values of the solid angle deficit, the 
effects induced by the conical structure and flux are exponentially suppressed. 



PACS numbers: ll.10.Kk, 04.62,+v, 98.80.Cq 



1 Introduction 

Symmetry braking phase transitions in the early universe have several cosmological consequences 
and provide an important link between particle physics and cosmology. In particular, within 
the framework of grand unified theories various types of topological defects are predicted to 
be formed pQ. Among them the cosmic strings are of special interest. Although the recent 
observational data on the cosmic microwave background radiation have excluded cosmic strings 
as seeds for structure formation, they are still candidates for the generation of a number of 
interesting physical effects such as the generation of gravitational waves and gamma ray bursts. 
Recently, cosmic strings attract a renewed interest partly because a variant of their formation 
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mechanism is proposed in the framework of brane inflation [2J. The non-trivial space-time around 
a cosmic string leads to interesting effects in quantum field theory. The vacuum polarization 
associated with scalar and fermionic fields in the geometry of an idealized cosmic string, have 
been analyzed in OH] and [5J, respectively. Moreover, considering the presence of the magnetic 
flux along the cosmic string, there appears an additional contribution to the vacuum polarization 
effect associated with charged fields [6]. The combined effects of the topology and boundaries 
in the geometry of a cosmic string are investigated in [TJ [8] . 

Though the topological defects have been first analyzed in four-dimensional space-time, they 
have been considered within the framework of the braneworld scenario as well. By this scenario 
our four-dimensional world emerges as a defect in a higher-dimensional space-time (for a review 
see [9]). Braneworlds naturally appear in the string/M-theory context and provide a novel setting 
for discussing phenomenological and cosmological issues related to extra dimensions. The models 
introduced by Randall and Sundrum are particularly attractive [10]. The corresponding space- 
time contains one or two Ricci-flat branes embedded on a five-dimensional anti-de Sitter bulk. 
More recently, the Randall-Sundrum scenario is generalized to the case of two extra dimensions 
by using a global string [llj . For the case with three extra dimensions, magnetic monopole and 
global monopole have been analyzed in [12] and [131 114j . respectively. In particular, in Q 1 the 
authors have obtained the solution to the Einstein equations considering a general p-dimensional 
Minkowski brane worldsheet and a d-dimensional global monopole, with d ^ 3, in the transverse 
extra dimensions with the core on the brane. 

The investigation of quantum effects in corresponding braneworld models is of considerable 
phenomenological interest, both in particle physics and cosmology. Quantum effects provide 
a natural alternative for the stabilization of the radion fields. The corresponding vacuum en- 
ergies give contribution to both the brane and bulk cosmological constant and, hence, has to 
be taken into account in the self-consistent formulation of the braneworld dynamics. In recent 
papers we have investigated the vacuum polarization effects associated with massless scalar [15] 
and fermionic [16] fields, respectively, in higher dimensional global monopole space-time, in 
braneworld context. One-loop quantum effects of a scalar field induced by a composite topo- 
logical defect consisting a cosmic string on a p-dimensional brane and a (m + l)-dimensional 
global monopole in the transverse extra dimensions are considered in [17] ■ In order to develop 
these analysis we have explicitly constructed the corresponding Green functions. The objective 
of this paper is to complete the analysis, studying the fermionic vacuum polarization effects 
associated with quantum charged field propagating on six-dimensional space-time produced by 
a composite defect, along the same line of investigation developed in [17] . Here we shall consider 
the monopole on the three-brane and the cosmic string on the two-dimensional transverse space. 
In this analysis we also allow the presence of an extra magnetic field, which can be understood 
as producing a magnetic flux running along the core of the cosmic string (about vacuum polar- 
ization energies for flux tube configurations see, for instance, [18j and references therein). Note 
that fluxes by gauge fields play an important role in higher-dimensional models including the 
braneworld scenario (see, for instance, [19] ). In particular, they provide an stabilization mech- 
anism for all moduli fields appearing in various string compactifications. The problem under 
consideration is also of separate interest as an example with gravitational and topology-induced 
polarizations of the vacuum, where all calculations can be performed in a closed form. 

This paper is organized as follows. In section [2j we introduce the structure of the space- 
time produced by the topological defect under consideration. With the objective to construct 
the fermionic propagator in this manifold, we first evaluate the heat kernel associated with the 
square of the corresponding Dirac operator. Although the general expression for the heat kernel 
is given in terms of series involving the modified Bessel functions, a much simpler expression 
is obtained for a particular choice of the parameters which codify the conical structure of the 
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transverse two-space and the fractional part of the ratio of the magnetic flux by the quantum 
one. This special case is considered in section [3j In this section we also analyze the fermionic 
propagator in the coincidence limit and explicitly extract the divergent part. By using this 
propagator, we evaluate the part in the vacuum expectation value of the energy-momentum 
tensor induced by the conical structure and the magnetic flux. In section H] we consider the 
general case for the parameters of the conical structure and the flux. Aiming to investigate 
the effects induced by this structure, we evaluate the subtracted heat kernel and the Green 
function. The corresponding part in the vacuum expectation value of the energy-momentum 
tensor is presented as a sum of two terms. The first one corresponding to the contribution 
coming from the bulk, having only a three-dimensional point-like global monopole on it, and the 
second one is induced by the conical structure of its two-dimensional submanifold. The behavior 
of the latter in various asymptotic regions of the parameters is discussed. The special case is 
discussed when the global monopole is absent. The main results of the paper are summarized 
in section [5l Throughout the paper the system of units h = c = 1 is used. 



2 Heat kernel 

The six-dimensional space-time that we want to investigate the fermionic vacuum polarization 
effects, is constituted by a point-like global monopole on a three-brane and a transverse conical 
two-dimensional space. In the coordinate system x A = (t, r, 9, (j), p, ip) it can be described by 
the line element 

ds 2 = g AB dx A dx B = -dt 2 + dr 2 /a 2 + r 2 (d6 2 + sin 2 9d<j) 2 ) + dp 2 + b 2 p 2 dip 2 , (1) 

with radial coordinates varying as r, p ^ 0, planar angles (f), <p € [0, 2ir], polar angle 9 € [0, tt] 
and t £ (— oo, oo). In (JT[) , the parameters a and b codify the presence of the global monopole 
and string respectively. In this space-time the scalar curvature outside the string axis is 1Z = 
2(1 - a 2 )/r 2 . 

The flat six-dimensional Dirac matrices, T^ M \ are 8x8 matrices, which can be constructed 
from the four-dimensional 4x4 ones |20| as shown below |21j : 

7 o")' r< M4 T)' r(5 ' = (-°; o- 

where p = 0, 1, 2, 3, 75 = Z7°7 1 7 2 7 3 , and / is the 4x4 identity matrix. These matrices obey the 
Clifford algebra {T^ M \ r^} = -2rf M W> I (8) , for M, N = 0, 1, ... 5. In this representation, 
the T^ 7 ) matrix takes a diagonal form below 

r (7) = r (o )r (D _ T ( 5) = n (3) 

and presents two chiral eigenstates defined as ^+ and [21J. So any six-dimensional fermionic 
wave-function, ^f, can be decomposed in terms of its chiral components as 

-(£)■ < 4 > 

Solutions for the Dirac equation, 

iT^d {M) * = M* , (5) 
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with defined chirality can only be possible for M = 0. In this way for positive chirality, equation 
dU) reduces to 

aWd {M) V + = 0, (6) 
being o^ M ^ = (7^, 175, —I) a set of 4 x 4 matrices; and for negative chirality it reduces to 



o^d iM) *^ = o, 



(7) 



being now <j( A/ ) = (7^, 175, I). 

For a massive charged fermionic filed in a curved space-time and in the presence of the 
electromagnetic field, the Dirac equation has the form 



(iJP + e4-M)V(x) = , 
with the covariant derivative operator defined by the relation 

W = r M (d M + n M ) , 

where 14^" is the spin connection given in terms of the T matrices by 



14 



M 



(8) 



(9) 



(10) 



and 



T M A 



M 



Now, in order to write the Dirac equation in the six-dimensional space-time defined by (JT]) 
we shall choose the following basis tetrad below: 
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For this basis tetrad, the only nonzero spin connections are: 



n e = -(l-a)f( g ) 

n * = "oC 1 - a) sin 0E (8 ) • 9. 



U„ = -(1-6)A 5 



(12) 



where 9 and (j) are the standard unit vectors along the angular directions on the three-brane, 



J (8) 



E 
£ 



with E 



Bp 
o P 



(13) 



with op being the Pauli matrices, and 

A - = i r (4) r (5) 



75 
-75 



(14) 
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The fermionic propagator defined by [22] 

iS F {x,x') = (0\T(^(x)^(x'))\0) , (15) 
with ^ = ^tr , satisfies the non-homogeneous linear differential equation, 

(iW + eA-M)S F (x,x') = ^=L6 6 (x-x') , (16) 

where g = det(gAB)- The Green function defined in (|16p is a bispinor, i.e., it transforms as ^ 
at x and as ^> at x'. 

Let a bispinor T> F (x,x') satisfy the differential equation 

[□ - ieg MN (D M A N ) + ieY, MN F MN - 2ieg MN A M V N 

-e 2 g MN A M A N — M 2 — 11/ A] V F (x,x') = -^=5 6 {x - x') , (17) 

with 

E MN = l[ r M )r JV] ; Dm = Vm _ ieAu } (lg) 

7£ being the scalar curvature and the generalized d'Alembertian operator given by 

□ = g MN V M V N = g MN (8 M V N + U M V N - { S MN }V S ) ■ (19) 

Then the spinor Feynman propagator can be written as 

S F (x, x') = {iJ7 + e4 + M) V F (x, x') . (20) 

Now we apply this formalism to the system under investigation. In order to take into account 
the presence of a magnetic field along the core of the cosmic string, we may write 

A M = Ad M y ■ (21) 
Choosing the basis tetrad (fTTj) . the operator /C, which acts on the bispinor in (fTT|) . reads 

K = -d 2 + a 2 (d 2 + 2 -d r y^ + d 2 + ^d p 

+ (d v + - ieA v f - (1 ~ a) [l + S (8) • lJ — M 2 , (22) 

where L is the ordinary angular momentum operator on the three-brane. As we can see, the 
above operator is expressed as diagonal blocks of 4 x 4 matrices. 

Because it is not possible to express the fermionic propagator for a massive field in the 
manifold given by ([I]) in terms of known special functions, we shall restrict our analysis for 
vacuum effects associated with massless field only. In this case we shall consider the field 
with positive chirality, for which the Dirac equation can be written in terms of a 4 x 4 matrix 
differential equation 

P*+ = , (23) 

with the operator 

P = 7°$ + \lr (ard r - S • L - 1 + a) + a ■ p (d p - ^) 



+^-<P {d v -ieA ), (24) 
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where a = (7-?, 275, —I)- In (|24[) we have r • a = f • 7 = 7 r , a ■ p = 275 cos 93 — Jsin</2 and 
a -(p = —275 sin — / cos being f , /S and <^ the unit vectors along the corresponding coordinate 
directions. The case of a field with negative chirality is considered in a similar way. 
The four-component Feynman propagator obeys the differential equation 

ipS F (x,x') = -1=5^ (x-x') , (25) 

and can be expressed in terms of the bispinor T> F by the relation 

S F (x,x) = ipV F (x,x') . (26) 
Now T> F (x,x') is a 4 x 4 bispinor, being the solution of the differential equation 

iCV F (x, x') = ^=5 (6) (x - x') , (27) 

with the operator 

JC = -d 2 + a 2 (d 2 r + 2 -d r ))-^ + d 2 + ±d p 

+ ^ + n <P ~ ieA f) 2 ~ ~2 1 + , (28) 

where 

tt v = 1(1 - 6)75- (29) 

The vacuum expectation value of the energy-momentum tensor can be expressed in terms of the 
Euclidean Green function, which is related with the ordinary Feynman Green function [22] by 
the formula D f (t, r; r', r') = —iT> F (x,x'), where t = it. In the following we shall consider the 
Euclidean Green function. 

In order to obtain the Euclidean Green function T> F (x,x') in explicit form, we need to have 
the complete set of normalized bispinors which obey the eigenvalue equation 

jC E ^\{x) = -X 2 ^x(x) with A 2 ^ , (30) 

being K E the Euclidean continuation of the differential operator given in (|28p . So, we may write 

V E (x, x') = £ ^ffifcO = r d§ £ * A(a0 *t (a /) . (31) 

The eigenfunctions of ([30]) can be specified by a set of quantum numbers associated with oper- 
ators that commute with iC E and among themselves. They are: p T = —id T , J 2 , J z , L 2 and S 2 , 
being S = and J = L + S, and p v = — id^ and 75. The latter two operators are associated 
with the variables on the conical two-space. 

On the basis of this set of operators, the normalized eigenfunctions of the operator K E are 
presented in the form 

A 2 = E 2 + a 2 p 2 + [3 2 , (33) 
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where represents the Bessel function of the orders 



j + l/2 n a 

v a = — , with n a = (-1) , cr = 0, 1 

a 2 

jig = n + 5^-^- - N - 7 , with 5 = +1. (34) 



In (|32p . tifj'/n- are ^ e s Pi nor spherical harmonics which are eigenfunctions of the operators I? 
and a • L as shown below: 



= + (35) 

with «(°) = -(Z + 1) = -(j + 1/2) and = Z = j + 1/2. In ([34} we have expressed 
e^4 = AT + 7, in terms of an integer number, TV, plus a fractional one, 7. As we shall see 
below, only the fractional part 7 leads to non-trivial effects. The index a specifies two types 
of eigenfunctions corresponding to I = j — n a /2 being I the orbital quantum number, and 
d = ±1 the corresponding eigenvalues of the 75 matrix. Moreover, we have E G (— 00, 00), 
n = 0, ±1, ±2, j = 1/2, 3/2, ... denoting the eigenvalue of the total angular quantum 
number, rrij = —j, j determines its projection and p, f3 G [0, 00). 
According to (|3ip the heat kernel is given by the expression 

/OO /*00 /'OO 

dE / dp dp V ^ 5)(CT) (x)^f (fT)t (x')e- sA2 . (36) 



Substituting (|32|) into the above definition, we get: 

1 „iiVA V / nn'\ 

K(x,x';s) = q , * K/0 YYliutUbl—i e inAlf 

V ; 32vr 3 / 2 6a V^ 7 s 5/2 V 2s 7 

where A 99 = ip — tp' , I^(z) is the modified Bessel function, 

q m .(o,Q0 = (38) 



is a 2 x 2 matrix and 



US = n + *(l-6)/2-7, 

V = a 2 (Ar 2 + p 2 + p /2 )+r 2 + r' 2 , (39) 



with At = t — t'. The general expression for the Green function is obtained by integrating ([3] 
as shown below: 

/•OO 

V E (x,x')= ds K(x,x';s). (40) 



For 6=1 and 7 = A" = 0, the summation over the quantum number n in ()3T[) can be done 
explicitly by using formula from [23] and the result coincides with the corresponding one given 

in ma. 



1 Explicit forms of above standard functions are given in [20] 
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3 Special case 



3.1 Green function 

Before to construct the Green function from the heat kernel (|37p for the general case of the 
parameters characterizing the conical structure and the magnetic flux, here we consider a very 
special case which allows us to obtain a much simpler expression. It has been shown in [3] that 
when the parameter q = 1/6 is an integer number, the scalar Green function in four-dimensional 
cosmic string space-time can be expressed as a sum of q images of the Minkowsiki space-time 
function. Also, recently the image method was used in [8] to provide closed expressions for the 
massive scalar Green functions in a higher-dimensional cosmic string space-time. The mathe- 
matical reason for the use of image method in these applications is because the order of the 
modified Bessel functions which appear in the heat kernel becomes an integer number. As we 
have seen, for the fermionic case the order of the Bessel function depends, besides on the integer 
angular quantum number n also on the the factor (1 — b)/2 from the spin connection. However, 
considering a charged fermionic field in the presence of a magnetic flux running along the string, 
an additional term will be present, the factor 7. In the special case 

7 = (1 - b)/2, (41) 

the order of the Bessel function becomes an integer number and the image method can be used 
to express the fermionic Green function in a closed form [21]. Here, the manifold that we are 
considering has the structure of a direct product of a cosmic string by a global monopole one. 
Accepting the above condition on the parameters b and 7, the calculation of the vacuum energy- 
momentum tensor becomes much easier to be performed as we shall see. So, although being a 
very special situation, the analysis of vacuum polarization effects in this circumstance may shed 
light on the qualitative behavior of these quantities for non-integer q. By using the formula [23], 

00 q-1 

I mg (x)e imqtp = Lj2e xcos(<p+27vk/q) , (42) 

m=— 00 ^ k=0 

after some intermediate steps, (l3?|) can be written as: 

K ^ X '^ = 327r3/2 a -f=, s 5/2 2-, 2^ M S 



rr k=0 " a=0,lj,mj 

(1 + e ^)q m .(o,Q') (1 - e^)ci m .(n,(V) 
(1 - e tf*)q; m .(fi,fi') (1 + e^)ci m .(n,n r ) 



(43) 



with notations 



14 = V -2a 2 pp' cos {Aip/q + 2kir/q) , 

P k = {l-l/q)A<p-2kn/q. (44) 

The corresponding Green function is obtained by substituting ([4"3]) into ([4*0]) . After the 
integration over the variable s with the help of formula from [23J, the result is: 

£l ' J 8vr 2 (r'r) 2 2^ 2^ 2^ sinhu 

^ ' fc=0cr=0,l j,mj K 



(1 - e^)C?m,n>) (1 + e^)Cf m 



(45) 
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with 

V fe _ a 2 (Ar 2 + p 2 + p' 2 - 2pp' cos(Aip/q + 2-nk/q)) + r' 2 + r 2 

COSU ti/^ — J — '. • 

2r r 2r r 

In (|45p . represents the associated Legendre function of the second kind. From the above 

expression we can see that, the k = component presents divergence at the coincidence limit; 
however, as to the others, they remain finite. The reason is because cosh Uk will be always greater 
than unity for these components, consequently the Legendre function assumes finite value. 



3.2 Vacuum average of the energy- momentum tensor 

The evaluation of the vacuum expectation value (VEV) of the energy-momentum tensor as- 
sociated with massless fermionic field in the six-dimensional space, in a braneworld scenario 
containing a three-dimensional global monopole, has been developed in [16]. Here we are mainly 
interested in the investigation of quantum effects induced by the presence of the cosmic string 
in the transverse two-dimensional space and by the magnetic flux along the core of the string. 

Using the point-splitting procedure [22], the VEV of the energy-momentum tensor associated 
with charged fermionic fields can be obtained by using the Feynmann propagator as shown below: 

{Tab(x)) = \ lim Tr [a A {D B - D B >) + a B {p A - D A ,)] S F (x, x') . (46) 

4 x — *x 

Here we shall use the four-component bispinor, where D B = Vb — iA B , the bar denotes complex 
conjugate, and a a = ej^ 0"(m)> with <j( M ) = (7^, ijs, I). Moreover, 

S F (x,x') =ipV F (x,x') , (47) 

being the operator ip given by (I24h . 

For the special case considered in this section, we shall use for the bispinor the result given in 
(|45p . As the first step in this direction, let us calculate the zero-zero component of (|46p . Because 
ao = — 7 , -Do = do and the dependence of the Feynman propagator on the time variable is by 
t — t' , we can write: 

(Too (»)) = lim Tr [ l0 d S F {x' ,x)) = - lim Tr [j d pV E (x', x)] , (48) 

x'— >x x'— >x 

where we have used the Euclidean version for the Green function. 

Before to embark in the above calculation it is useful to analyze, separately, the contribution 
coming from the k = component of the Euclidean Green function. As the first point we can 
see that taking 0, = 0' in CJ m , the summation over the quantum number rrij in (1451) provides: 

Y. C U^) = 2J ^h) ■ (49) 

mj 

On the other hand, taking ip' = ip, the k = component of the Green function becomes 
proportional to I( 4 ). The same is true when the operator S • L of (f24"|) acts on CJ m ., because of 
(|35p . In this way, almost all derivative terms provide a vanishing trace. The exception is for the 
term proportional to cL; however, after taking this azimuthal angle derivative on (|45p and the 
coincidence limit ip' = 93, the matrix left is proportional to i(l — b){I^ — 75), consequently its 
products with 7570 and 70 have a vanishing trace. As a conclusion we can affirm that the only 
term that produces a non- vanishing result in the calculation of the vacuum average is linear in 
the time derivative in (|24|) . consequently for the k = component we have: 

<T 00 (*)> (fc=0) = - lim Tr d 2 V^ =0 \x, x') . (50) 
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Taking in T>g~°\x, x') the coincidence limit for the azimuthal angle in the extra conical two- 
dimensional space, this component coincides with the corresponding bispinor in the absence of 
string and magnetic flux. As a consequence the analysis of the above vacuum average is similar 
to that one developed in [16J. The new analysis will involve the k ^ 1 components of (|45j) . As 
we have already mentioned, the main objective of this paper is to evaluate the quantum effects 
induced by the string and magnetic flux on the VEV of the fermionic energy-momentum tensor. 
So, our next step is to calculate the contributions from the string and magnetic flux to (Tqq(x)). 
We will denote this part as (Tqq(x)) s . 

Because the k 1 components of (|45p are finite at the coincidence limit for p ^ 0, the 
calculation of their contribution to the VEV of the zero-zero component of the energy- momentum 
tensor does not require any renormalization procedure. Moreover, this contribution is only due 
to the term of the operator (|24p with time derivative . The reason resides in the fact that only 
second time derivatives of 1/ sinhufc and (^^^(cosh u&) in (135]) . present non- vanishing result 
in the coincidence limit. So for these component of the bispinor we also can write: 

(r 00 (x)) s = - Km T* [^dopV^ix,^ = - Jim Tr d 2 T V { * >l \x,x') . (51) 

Now, after some intermediate steps and introducing the rescaled radial coordinate r = r/a, we 
arrive at the expression: 

^°>s = -^E t^^ D [Ql/a-li^+Ql/M] , (52) 
Ibir^r^y f— * 1 + y sin {TTK/q) ~: L ' 'J 

where 

y = (p/r) 2 , z k = 1 + 2ysm 2 {'Kk/q) . (53) 

By making use of the recurrence relation for the associated Legendre function of the second kind 
(see [25J), it can be seen that 



Ql-i(z) + Ql(z) = -2yJ [(1 + v)Ql_ x {z) + (1 - v)Ql(z)] , (54) 

and, hence, the summand in (j52H can also be written in terms of the associated Legendre function 
of the order 1. 

When the global monopole is absent one has a = 1, and by using the formula 

oo 

J2(2l + l)Q?(x)=2j^—2, (55) 
i=o \ ' 

from (1521) we find 



/t0\ 1 ^ cos 2 (7rfc/g) q 2 - 1 2 2 

(To)s = -lW g sin«(W?) = "^rV (9 " 4)((? +5) - (56) 

In section |3] we shall see that these result holds for non-integer values of q as well. By taking 
into account that the same result will be obtained for the negative chirality, in the case of the 
combination of both chiralities formula (|56p coincides with the result given in [23] . 

For q = 2, the expression in (I52D vanishes, and there is no contribution to the VEV of the 
energy density induced by the cosmic string and magnetic flux. However, for q = 3 there appears 
a non-vanishing contribution given below: 

—2 i 00 

(T ° 0)s = - 2A^yl + 3y/4 E ' [<#/«-! ^ + 3 ^ 2 ) + C 1 + ^ " ^ 57 ) 
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In figure [Q we have plotted the behavior of f 6 (TQ) s , evaluated by formula (|57p. as a function 
of the parameters a and y. As it is seen from this figure, the part in the VEV of the energy 
density tends to zero for small values of the parameter a. The behavior of the VEV in various 
asymptotic regions of the parameters will be described below for the general case when there is 
no relation between the planar angle deficit and the fractional part of the magnetic flux. 




Figure 1: The energy density induced by the string, f 6 (TQ ) s , as a function of a and y = (p/r) 2 
for q = 3, 7 = 1/3. 

In figure [21 we have presented the dependence of the part in the energy density induced by 
the string, on y = (p/r) 2 for different values of the parameter q (the numbers near the curves). 
The corresponding values of the parameter 7 are found from condition (I41jl . The graphs are 
plotted for the value of the global monopole parameter a = 0.5. Note that in the cases q = 1, 2 
the corresponding energy density vanishes. 

The VEV of the radial component of the energy-momentum tensor will be given by 

(T rr (x)) = -!- lim Tr [f ■ 7 (d r - fy)] S F (x, x') , (58) 
la x'-^x 

where we have used D r = d r and a r = r ■ j/a. Expressing the Feynman propagator in terms of 
the Euclidean Green function, we can write: 

(T rr (x)) = -— lim Tr [f ■ j (d r - d r >) pV E (x\ x)\ . (59) 
la x'^x L J 

At this point we can make a similar analysis as we did in the calculation of the VEV of the 
energy density. For the k = component of (145 j) . the only non- vanishing contributions coming 
from the product of 7 r by (|24|) . are given by the second, third and fourth terms of the latter. So, 
we may take for this component <p' = <p. As a consequence the corresponding vacuum average 
again coincides with the one in the absence of string and magnetic flux. So, let us now proceed 
the calculation of the new contributions for the above VEV due to the terms with k ^ 1 in (1451): 



{T rr (x)) s = --!- lim Tr \r -j(d r - d r >) pV% >K \x,x')] . (60) 
la x'—*x L 
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0.2 0.4 0.6 0.8 1.0 

(p/r) 2 



Figure 2: The dependence of the energy density induced by the string, f^{T^) s /q, as a function 
of y = (p/r) 2 for a = 0.5. The numbers near the curves correspond to the values of the 
parameter q. 



Because of some simplifications that appear in this calculation, it is possible to observe that 
the only non-vanishing contribution comes from the second term on the right-hand side of (j2 
and one finds 



1 



{T rr (x)) s = - lim Tr 



d r (d r — d r i) 



a 



2 x'^x 

Now substituting ([4"5]) into (|6ip , and using relation ([33]) , we can see that 

(T r r ) s = (T °) s . 

The calculation of the component (T pp ) is formally given by 



1 



ppi 



- lim Tr [a p (d p - d p >)] S F (x,x') , 



(61) 



(62) 



(63) 



where we have used D p = d p . Expressing the Feynman propagator in terms of the Green 
function, and analyzing the k = component of (|45p . we conclude, after a careful analysis, that 
this component provides the same vacuum average as in the absence of string and magnetic flux. 
The only non-vanishing contributions to the parts in the VEV due to the string and magnetic 
flux come from the last two terms on the right-hand side of (|24p : 



1 



{TP(x)) s = - lim Tr 



d p (d p - d pl ) + L_L _ ™.(id v + JV + 7 ) 



Pg^(x,x'), (64) 



where 7 = (q — l)/2q. Substituting the expression for the Green function, and after a long 
calculation, we obtain for VEV of T p the same expression as obtained for Tq, given in f)52[) : 



(TPs = (T °) f 



(65) 



The other components of the VEV for the energy-momentum tensor are found from the 
covariant continuity equation V^(Tg) = 0. For the geometry under consideration this equation 
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reduces to the following two relations 

(Tt) = (T+) = (l + r -d r ) (2*) , 

(T$) = 3>(r;)), (66) 

By using formulae (j62[) . ()65|) . (I66h . it can be checked that the part in the VEV of the energy- 
momentum tensor induced by the string and magnetic flux is traceless, (T^) s = 0, and the trace 
anomaly is contained in the pure global monopole part only. 

Having obtained the complete expressions for the VEVs of all components of the fermionic 
energy-momentum tensor induced by the cosmic string and magnetic flux, we let for the next 
section the analysis of the general case, where there is no more relation between the magnetic 
flux and the parameter associated with the cosmic string. 



4 General case 
4.1 Green function 

In order to evaluate the vacuum expectation value of the energy-momentum tensor induced 
by the magnetic flux and cosmic string for general values of the parameters b and 7, we shall 
subtract from the heat kernel (|37p the corresponding heat kernel in the absence of cosmic string 
and taking 7 = 0. H Denoting by /C au b(x, x'\ s) the subtracted heat kernel and introducing 
rescaled coordinates (p = b(p and r = r/a, we have 



327r 3/2^~~7 s 5/2 

<t=U,1 J,rrij 



*EE 



J,m,j \ 1 ) 3, Taj 




(67) 



where Aip = (p — dp' and 

V = Ar 2 + p 2 + p' 2 + f 2 + r 2 . (68) 

Though the contribution of the separate terms in square brackets of (I67D to the Green function is 
divergent in the coincidence limit, the part of the Green function corresponding to the subtracted 
heat kernel (j67[) is finite in the coincidence limit for points away from the string core. This follows 
also from general arguments. For points outside the string core the electromagnetic field is zero 
and the local geometry is the same as that in the case when the global monopole is present only. 
Hence, the divergences in the Green function in the coincidence limit are the same in these two 
situations. 

To provide a more convenient expression for the subtracted heat kernel, for the summation 
over n we apply the Abel-Plana formula in the form 

E F(n±P)= r d uF(u) + i i^du £ gj^g^p (69) 

n=n± J0 J0 A=±l 



2 It is possible to see that the expression obtained from (|37[) by taking b — 1 and 7 = coincides, up to the 
factor e tNA f with the expression for the heat kernel given in [16] 
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where n + = 0, n_ = 1, and < (5 < 1. This generalization of the Abel-Plana formula with the 
upper sign is given in [26J and the formula with the lower sign is easily obtained from the upper 
sign case by redefining n — > n + 1. Note that formula (|69|) is also obtained from the generalized 
Abel-Plana formula (see |27J). 
Defining the sum 



JNqAp 



5> 



q\ns\ 



z e 



iqnAtp 



(70) 



and using formula (|69p . after some intermediate steps we find 



2e 



i(N-/3 s )qA<p 



1 



H — / dusmh(u7r)Ki u (z) 

J v \ L 



dul u (z) cos{uA(p) 

e -sgn(l3 s )\uA(p 



A=±l 

where K v (z) is the MacDonald function and we have defined 



o2ir{u/q+i\\l3 6 \) _ \ 



(3 S = S(q - l)/2q - 7 . 



(71) 



(72) 



To the sum over n with the second term in the square brackets of (|67p we apply the Abel-Plana 
formula in the standard form. This gives 



S(z,A<p) 



Em 



z e 



inA(p 



7T 



H — / dusmh(iru)Ki u (z) 



dul u (z) cos(uAif) 
cosh(uA(f) 



1 



(73) 



The subtracted Green function, 2? su b(x, x'), that we shall use to calculate the vacuum polar 
ization effects induced by the cosmic string and magnetic flux, is obtained from formulae (|40l 
and (E 



V suh (x,x') 



V2a~ 2 



167TV2, 



rr Jo 



du « 1 / 2 e-«V/2 



^Sj$(u P f/,A<p)-S(upf/,A<p) 



cr=0,l j,mj 



q w .(o,o') <sc£ m .(fi,n') 



(74) 



where V is defined by relation (|68p . The divergent contributions from the separate terms in the 
square brackets of (|74p come from the first integrals on the right-hand sides of formulae (|7ip 
and (173p . Now we see that after the application of the Abel-Plana formula these contributions 
are explicitly cancelled out in the subtracted Green function for A(p = 0. 



4.2 Energy-momentum tensor 

Having the subtracted Green function we can evaluate the VEV of the energy-momentum tensor 
on the base of formula (|46p . As in the special case discussed before, let us consider the part in 
the vacuum average of Tg induced by the string and flux: 

(T °(x)) s = lim Tr [ lQ d pV auh {x,x')\ . (75) 

X— >x 

Being P su b(a;, x') finite at the coincidence limit, we can interchange the differential operator and 
the integral over u. Moreover, terms linear in time derivative acting on the exponential factor 
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of the subtracted heat kernel, e _V//4s , produce a term linear in At, which goes to zero at the 
coincidence limit. The only term that survives is the one obtained by the second time derivative. 
So on the basis of these information, we may write: 



{T°(x)) s = hm Trd 2 T V suh (x,x') = lim Tr / ds d 2 T 1C snh {x, x'; s) 

x— >x x'^x Jq 



(76) 



Developing all the calculations needed and defining a new integration variable, v = uf 2 , we 
obtain: 

(7=0,1 j 

x^WE ^g/,| w |W-24W • (77) 

n <-g=±l ■ 

with y defined by ([551) . 

An equivalent form for (Tg ) a is obtained by using formulae (|7ip and (|73p with Aip = 0: 



/n />00 00 

< T o°>s = - l6 J /2a 2,e I dv v*' 2 e-^y> ^ I [ + W/ 2 *^M, (78) 



x/2 

— i avv-e — 

i=l 

In this formula we have introduced the notation 



^OO = J2 SjfUz,0)-2S(z,0) = - / dusmh(uTr)K iu (z)G q>J (u). (79) 
i=±i ' * Jo 

where ^ 

Gq n {u) = ^ ^ p 2Tr(u/q+i\\P s \) _ 1 ~~ p 27ru _ 1 ' ( 80 ) 

5=±l A=±l 



and /?5 is defined by (|72p . Note that, as it follows from (|80p . the VEV of the energy- momentum 
tensor does not depend on the sign of the fractional part of the magnetic flux. 

We could have obtained the second relation in ([79]) by applying to the series over n in the 
definitions of Sg\j (z, 0) and S (z, 0) the summation formula 



00 i-oo roo 

V F(\n + P\) = 2 / duF(u)+i V 



F{iu) - F(-iu) 



A=±l C 



which directly follows from (|69p . Note that for this formula the restriction < j3 < 1 is not 
necessary and it holds for any real value of the parameter (3. The Abel-Plana formula in its 
standard form (see, for instance, [28]) is obtained from here taking (3 = 0. 

The components (T rr ) s and (T pp ) s are found by the formulae similar to ([6"Tj) and ([M]) 

with the replacement T>^^ l \x,x') — ► 2? su b(x, re'). The resulting expressions for (T,T) S and 
(Tp) s are obtained from formula flTTJ for the energy density by the replacement I Vg (z) — ► 
[Iy a +n a {z) + I Ua (z)]/2 with z = f 2 /2s in the case of (T^) s and by the replacement Igi^i (z) — > 
[Iq\iis\+6p s /\ps\ ( z ) ( z )]/2 with z = p 2 /2s in the case of (Tp) s . After summation over <r and 

we see that (T r r ) s = (T p p ) s = (T °) s . Other components are found from relations (166p . It can 
be additionally checked that the resulting tensor is traceless. 

Now we consider the behavior of the string induced part in the VEV of the energy density in 
the asymptotic regions of the parameter y. For large values of this parameter, y>l, introducing 
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in (|78[) a new integration variable z = yv and expanding the integrand over 1/y, to the leading 
order we find 

(i ° )s ~ 42+1/^3^^6 1 duG ^^r(l/a + l/2)r(l/a + 5/2)' (82) 

for /5>f. In particular, for a < 1 the expectation value induced by the string and magnetic flux 
vanishes on the core of the global monopole. From formula (|82p it follows that at large distances 
from the string, the string induced energy density is suppressed by the factor (r/ p) 2 l a ~ 2 with 
respect to the corresponding quantity in the absence of the global monopole. 

As the vacuum expectation value given by (|78p diverges on the string core corresponding to 
y = 0, in the limit y <C 1 the main contribution into the sum over I comes from large values 
I and we can use the uniform asymptotic expansion for the functions Ii/ a ±i/2{'v). As the next 
step, we replace the summation over I by the integration. After some intermediate calculations, 
for the vacuum average of the energy density induced by the presence of the string and magnetic 
flux to the leading order one finds 

1 r°° 

^ w -RfTTT / duu{u 2 + l){u 2 + A)G^(u). (83) 

The expression on the right-hand side of this formula does not depend on the parameter a and 
coincides with the corresponding quantity in the geometry of a cosmic string with magnetic flux 
when the global monopole is absent (a = 1, see below). This means that near the cosmic string 
the most relevant contribution to the vacuum expectation value comes from the string itself. 

For small values of the parameter a, corresponding to strong gravitational field, the order 
of the modified Bessel functions in (|78p is large. Again, replacing these functions by the cor- 
responding uniform asymptotic expansion we can estimate the integral over v by the Laplace 
method. In this way it can be seen that the main contribution comes from I = 1 term and one 
finds 

„ exp[-(l/q)ln(l + 2y + 2^(l + y))] f°° 

As we see, in this limit the vacuum expectation values are exponentially suppressed. The similar 
feature takes place for the VEV of the energy-momentum tensor and the fermionic condensate 
induced in the global monopole bulk by the presence of boundaries on which the fermionic field 
obeys the MIT bag boundary condition [29J. 

In the case when the global monopole is absent, a = 1, the summation over / in the expres- 
sion ([75]) for (Tq) s can be done by using the formula [23] ^2^ Q (2l + l)^/+i/2( x ) = e x \f2x/TT. 
Redefining a new variable z = yv we have 

1 f'°° 

(T °° )s = "8^7 J Q dz z 2 ^^) ■ (85) 

By taking into account the expression (|79p for F q ^(z), the integral over z is evaluated explicitly 
using the formula from [23] . After some intermediate steps, we arrive to the following expression: 

1 

(r °° )s = "60^6 J o duu{u 2 + l){u 2 + A)G^(u). (86) 
Note that by making use of the formula 

poo r | 

A ^P) = / dun 2 -' E^R^ 

•^0 \ LI 



(u/q+iXf3) _ 1 p2iru _ ]^ 

,A=±1 ~ 



(_-t\m~l 
V L ) r 2m 



2m 



[q 2m B 2m ((3) - B 2m ] , (87) 
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where i?2 m (/3) an d ^2m are the Bernoulli polynomials and Bernoulli numbers respectively, the 
energy density is also presented in the form 



In the presence of both chiralities the result given by (|89|) should be doubled and it coincides 
with the formula derived in [23] . Considering the special case given by (|4ip , from (|89p we obtain 
formula (|56p . Hence, though we have derived formula (|56|) for integer values of q, we see that it 
is valid for non-integer q as well. 

5 Conclusions 

In this paper we investigate quantum vacuum effects for a massless fermionic field induced by a 
composite topological defect in a six dimensional spacetime. The spacetime is a direct product 
of the two dimensional cosmic string and three dimensional global monopole geometries. The 
corresponding heat kernel is constructed and on the base of this, the Green function is evaluated 
for a field with positive chirality. The consideration of a negative chirality field is similar. As 
the corresponding geometry with a global monopole and in the absence of a cosmic string was 
investigated previously in [16j . here we mainly concentrate on the effects induced by the cosmic 
string and magnetic flux along the core of the string. With this aim, the complete Green 
function is presented as the sum of two terms. The first one contains information only on the 
global monopole defect and the second one is induced by the presence of the cosmic string and 
magnetic flux. For points away from the string the second term in the Green function is finite 
in the coincidence limit and can be directly used for the evaluation of the vacuum expectation 
value of the energy-momentum tensor. 

In section [3] we have considered a special case when the parameters of the cosmic string and 
magnetic flux are connected by relation (|4ip and q is an integer. In this case the summation 
over n in the mode-sum for the heat kernel can be done explicitly and the Euclidean Green 
function is presented in the form (145 p . The vacuum expectation value of the energy- momentum 
tensor is constructed from the Feynman propagator by making use of formula (|46f) . The energy 
density induced by the string and magnetic flux is given by formula (|52p . The corresponding 
vacuum stresses along r- and p-directions coincide with energy density. The other components 
of the vacuum stress are found from the covariant continuity equation for the energy-momentum 
tensor. We have explicitly checked that the part of the energy- momentum tensor induced by 
the string and flux is traceless. 

The Green function and the vacuum expectation value of the energy-momentum tensor for 
the general case of the string and flux parameters are discussed in section SJ By using the Abel- 
Plana summation formula, the corresponding subtracted Green function is presented in the form 
()74p where the parts giving the divergences in the coincidence limit are explicitly cancelled out. 
Further, this Green function is used for the evaluation of the corresponding energy-momentum 
tensor. The vacuum energy density is given by formula ()78p and, as for the special case, the 
stresses along r- and /3-directions coincide with energy density. The remained components are 
found from relations (j66H . As for pure cosmic string geometry, the VEV of the energy- momentum 






(89) 
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tensor depends only on the fractional part of the magnetic flux and is an even function of this 
parameter. We have investigated the vacuum energy- momentum tensor in asymptotic regions 
of the parameters. At large distances from the string, p 3> f , the leading term of the asymptotic 
expansion of the energy density is given by formula (|82p . In this case the string induced energy 
density is suppressed by the factor (r/p) 2j/a ~ 2 with respect to the corresponding quantity in the 
absence of the global monopole. For points near the string core, p <C f, to the leading order 
the string induced part in the vacuum average coincides with the corresponding quantity in the 
geometry of a cosmic string with magnetic flux when the global monopole is absent. For small 
values of the parameter a, corresponding to strong gravitational fields, the vacuum expectation 
values are exponentially suppressed. 
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